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Abstract. We calculated a grid of evolutionary models for 
white dwarfs with helium cores (He-WDs) and investigated 
the occurrence of hydrogen-shell flashes due to unstable hy- 
drogen burning via CNO cycling. Our calculations show that 
such thermal instabilities are restricted to a certain mass range 
(M w 0.21 . . . 0.30 M0), consistent with earlier studies. Mod- 
els within this mass range undergo the more hydrogen shell 
flashes the less massive they are. This is caused by the strong 
dependence of the envelope mass on the white dwarf core mass. 
The maximum luminosities from hydrogen burning during the 
flashes are of the order of 10^ Lq. Because of the develop- 
ment of a pulse-driven convection zone whose upper boundary 
temporarily reaches the surface layers, the envelope's hydrogen 
content decreases by AAT « 0.06 per flash. 

Our study further shows that an additional high mass-loss 
episode during a flash-driven Roche lobe overflow to the white 
dwarf's companion does not affect the final cooling behaviour 
of the models. Independent of hydrogen shell flashes the evo- 
lution along the final white dwarf cooling branch is determined 
by hydrogen burning via pp-reactions down to effective tem- 
peratures as low as « 8000 K. 



are based on a 1 Mq model sequence extending from the pre- 
main sequence stage through the red-giant branch (RGB) do- 
main. We forced the models to move off the giant branch and to 
evolve into the white dwarf regime by applying large mass-loss 
rates at appropriate positions to take into account the binary na- 
ture of He-WDs (for details see Paper I). 

As pointed out in Paper I one of the major results of our 
study was the dominant contribution of hydrogen burning to 
the luminosity budget of the He-WDs. Therefore the final cool- 
ing evolution is slowed down, and the derived cooling ages are 
notably larger than those found in models which do not con- 
sider nuclear burning or do not find hydrogen burning to be 
important due to much lower envelope masses. 

In the present paper we will discuss the evolution of se- 
quences which undergo hydrogen shell flashes in detail. It is 
organized as follows: In Sect. 2.1 we will briefly summarize the 
main reasons for unstable nuclear burning, and report in Sect. 

on former studies concerning instabilities in He-WDs. In 
Sect. H we describe the evolutionary code used for our calcu- 
lations. The main results of the calculations are discussed in 
Sects. [4-.l| to 4.4 . Finally, conclusions are given in Sect. |. 
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1. Introduction 

In Driebe et al. (1 998) (hereafter referred to as Paper I) we pre- 
sented a grid of evolutionary tracks for low-mass white dwarfs 
with helium cores (He-WDs) in the mass range from 0.179 to 
0.414 Mq. The lower masses allow applications to compan- 
ions of millisecond pulsars. As an example we derived a cool- 
ing age for the He- WD companion of the millisecond pulsar 
PSR J1012H-5307 of ^ 6 Gyr, which is in good agreement with 
the pulsar's spin-down age of ~ 7 Gyr. The evolutionary tracks 
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2. Unstable nuclear burning 

2.1. The two general cases 

There are two reasons for nuclear burning to become unstable: 
The first one is the decoupling of the thermal and mechanical 
structure of a star due to large electron degeneracy. If the equa- 
tion of state gets more and more independent of the tempera- 
ture a local increase in energy production cannot be stabilized 
by a local expansion with following cooling of the affected lay- 
ers. Therefore a further increase in energy production will lead 
to a corresponding increase in temperature which in turn will 
again raise the energy release. This phase of thermally unstable 
burning will continue until the ongoing temperature increase 
results in an effective lifting of degeneracy in the burning zone. 
The stronger coupling of thermal and mechanical structure will 
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then allow for an expansion and the transition to a phase of 
stable nuclear burning. This kind of instability is found, for in- 
stance, during the onset of central helium burning in low-mass 
stars (A/initial ^ 2M0, A/coro « 0.47 Mq for Z = 0.02) on 
the tip of the red giant branch (central helium flash). 

The second reason for thermally unstable nuclear burning 
is mostly based on the geometry of the burning region and oc- 
curs only during shell burning. If the shell's mass (its "thick- 
ness") becomes too small compared to its radial extent, the ex- 
pansion following a local increase in temperature and energy 
production is insufficient to cool the shell: The thin burning 
zone will instead be heated by expansion, and the energy pro- 
duction is further increased. The thermal runaway only stops 
when the thickness of the shell is large enough as to allow for 
a cooling expansion. This kind of unstable burning repeatedly 
takes place during the double shell-burning phase of AGB stars 
(Schwarzschild & Harm 1965, Weigert 1966). There, these in- 
stabilities are known as thermal pulses or helium shell flashes 
caused by thermally unstable helium burning. 

Schwarzschild & Harm (1965) derived an instability crite- 
rion for non-degenerate matter by linear perturbation analysis. 
This criterion includes the shell thickness as well as the tem- 
perature exponent v of the energy generation rate ( e oc T^) as 
a measure of the temperature dependence of nuclear burning. 
The assumption of non-degeneracy is justified for the helium 
layers of AGB models during most of the thermal pulse evolu- 
tion. According to the Schwarzschild & Harm criterion unsta- 
ble burning is favored when the shell becomes thinner and the 
burning is quite temperature sensitive, as it is the case for he- 
lium burning. A similar criterion for the study of thermal pulses 
was derived by Sackmann (1977). In a recent study Frost et al. 
( 1 998) report that in advanced stages of thermal pulse evolution 
degeneracy in the region of the helium burning shell may be- 
come noticeable and can therefore lead to significantly stronger 
pulses (degenerate pulses). Kippenhahn & Weigert (1990) give 
an instability criterion which also takes degeneracy into ac- 
count. 

2.2. Unstable burning in helium white dwarfs 

Unlike as in AGB stars, where the helium burning shell be- 
comes thermally unstable, He-WDs show instabilities related 
to CNO cycling which dominates hydrogen burning in the 
lower, i.e. hotter part of the geometrically thin shell. 

Kippenhahn et al. (1967) calculated a He- WD model with 
M — 0.264 M0 which evolves through a phase of unstable 
burning. As a result the track in the Hertzsprung-Russell dia- 
gram (HRD) is reversed and the He- WD returns to the RGB 
domain. The evolution of the white dwarf during the thermal 
instability was followed in more detail by Kippenhahn et al. 
(1968). Giannone et al. (1970) found unstable burning in a 
0.268 Mq He- WD model, whereas their models with M — 
0.366 Mq and 0.426 Mq did not show any sign of thermal in- 
stabilities. Unlike the previously mentioned studies, Webbink 
(1975) found only small hook-like excursions in his He- WD 
tracks close to the point of maximum effective temperature. 



Iben & Tutukov (1986) calculated a 0.298 Mq He- WD 
model based on a 1 Mq model which suffered from high mass 
loss episodes on the RGB, mimicking Roche-lobe overflow to a 
companion. They found two strong thermal instabilities to oc- 
cur on the cooling branch resulting in a track similar to the one 
of Kippenhahn et al. (1968). 

Castellani et al. (1994) calculated a grid of He-WD model 
sequences for different metallicities Z using the same tech- 
nique as Iben & Tutukov (1986). Like Giannone et al. (1970) 
they found thermal instabilities only below a certain mass limit 
wich depends on Z: For Z = 0.01 and 0.001 no flashes 
were found for M k, 0.33 Mq and M ^ 0.35 Mq, respec- 
tively. For Z = 2 ■ lO""* one hydrogen shell flash was found 
for M = 0.370 Mq and M = 0.389 M©. In a recent study 
Sarna et al. (1998) found instabilities in their cooling tracks for 
M < 0.2 Mq, quite similar to those of Webbink (1975) who, 
however, used rather large time steps. 

Due to large time steps and/or significantly smaller enve- 
lope masses several studies did not find thermal instabilities at 
all during the cooling of He-WDs, as e. g. Chin & Stothers 
(1971), Alberts et al. (1996), Althaus & Benvenuto (1997), 
Benvenuto & Althaus (1998) and Hansen & Phinney (1998). 

3. The evolutionary calculations 

Besides some minor modifications we used the evolutionary 
code described by Blocker (1995). Nuclear burning is ac- 
counted for via a nucleosynthesis network including 30 iso- 
topes with aU important reactions up to carbon burning sim- 
ilar as in El Bid (1994). The most recent radi ative opacities 
by Iglesias et al. (1992) and Iglesias & Rogers ( 1996 ), supple- 
mented by those of Alexander & Ferguson (1994) in the low- 
temperature region, are employed. Diffusion is not considered. 
The initial composition is {Y, Z) ~ (0.28,0.02), the mixing 
length parameter a = 1.7 followed from calibrating a solar 
model. The Coulomb corrections to the equation of state are 
those given by Slattery et al. (1982). 

He-WDs are known to be components in binary systems 
where early case B mass transfer (Kippenhahn & Weigert 
1967) must have taken place during RGB evolution which 
forces the stars to leave the RGB before the onset of helium 
burning and to evolve into the white dwarf regime. We did 
not calculate the mass exchange phases during the RGB evolu- 
tion in detail because we are primarily interested in the cooling 
properties of the white dwarf models themselves rather than in 
the generation of these models by binary star evolution. Hence 
we used an approximate approach to get the pre-white dwarf 
models (see also Iben & Tutukov 1986, Castellani et al. 1994, 
and Paper I): 

We calculated a 1 Mq sequence from the pre-main se- 
quence phase up to the tip of the RGB. Along the RGB we 
applied mass-loss rates according to Reimers (1975) with 
rj — 0.5. At appropriate positions high mass loss rates, Afhigh, 
were invoked in order to get models of desired final mass, M: 
0.179, 0.195, 0.234, 0.259, 0.300, 0.331 and 0.414 M©. Fig- 
ure |l] displays aU He-WD sequences in the HRD. High mass 
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Fig. 1. Hertzsprung-Russell diagram with evolutionary tracks 
for He-WD models. The tracks with solid lines belong to final 
white dwarf masses of 0.179, 0.195, 0.234, 0.259, 0.300, 0.331 
and 0.414 Mq (from top to bottom). The dotted lines show the 
flash phases for M = 0.234 and 0.259 M©, the dashed fine 
gives the IMq evolutionary path the He-WD models are based 
upon. 

loss varied from Mhigh ~ 10"^ Mq yr^^ for M w 0.15 Mq 
to about 10~^ Mq yr^^ for M « 0.4 M©. These values were 
chosen to allow the models to hold thermal equilibrium during 
their further evolution (with M — Afhigh)- For a more detailed 
description we refer to Paper I. Here, we only want to stress 
the importance of sufficiently small time steps when the evolu- 
tion through thermal instabilities should be followed properly. 
A brief discussion on this topic is given in the Appendix. 

4. Results 

4.1. General remarks 

We briefly repeat the main results of our calculations which 
have been addressed in detail in Paper I. Due to the compar- 
atively large remaining envelope masses after termination of 
the RGB evolution (see Table [l|]) one of the main character- 
istics of our white dwarf models is that hydrogen burning due 
to pp-reactions remains the main energy source down to effec- 

' We note that the same table is shown in Paper I but the data for 
M = 0.179 and 0.195 Mq did not refer to T^s = 10000 K as indi- 
cated in the caption but erroneously to Toff — 5000 K as for the larger 
masses. 



Table 1. Total remnant mass M, mass of the hydrogen- 
exhausted core Mc, total mass of the outer hydrogen layers 
("thickness") A/h, envelope mass Alcnv, and helium surface 
abundance by mass fraction Y, at Tcs = 5000 K for M > 
0.2Mo and at 10000 K for M < O.2M0 after flie end of RGB 
evolution. Mc is defined by the mass coordinate below which 
X < 0.35 with X being the mass fraction of hydrogen. 



M/Mq 


Mc/Mq 




Afonv 

IO-^'Mq 


Y 


0.179 


0.1693 


5.061 


10.211 


0.464 


0.195 


0.1859 


4.937 


9.598 


0.461 


0.234 


0.2220 


8.118 


13.098 


0.354 


0.259 


0.2524 


4.771 


7.232 


0.312 


0.300 


0.2960 


3.189 


4.746 


0.301 


0.331 


0.3281 


2.509 


3.744 


0.301 


0414 


04116 


1.446 


2.175 


0.301 



tive temperatures well below lO^'K (see Paper I and Fig. ||). 
This residual burning leads to a significant slow-down of the 
further evolution, resulting in larger cooling ages (typically a 
few 10^ yr, see Fig. ^) than found for He-WD models where 
hydrogen burning is negligible due to smaller, ad-hoc assumed 
^ envelope masses, or even not considered at all. The impli- 
cations of evolutionary envelope masses for the evolution of 
white dwarfs are discussed in Blocker et al. (1997). 
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Fig. 2. Ratio of hydrogen (inyd) to gravothermal luminosity 
(ig) as a function of Tcff for He-WDs of different masses (in 
Mq). For the sequences which undergo hydrogen flashes only 
the final evolution along the cooling branch is shown. 

Besides of this important property of He-WDs and its con- 
sequence for age determinations of, e.g., millisecond pulsar 
systems as PSR J1012H-5307 (cf. Paper I), another result of 
our study concerns the occurrence of hydrogen shell flashes: 

^ In these calculations the envelope mass is taken as a free parame- 
ter and has not been computed from the mass-loss history of the He- 
WD progenitor according to the binary nature of He-WDs. 
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by the fading CNO cycling (see next section) is restricted to a 
certain mass range of M « 0.21 . . . 0.30 M0. Furthermore, the 
number of flashes increases with decreasing white dwarf mass. 

It is noteworthy that the cooling properties below Teff ^ 
20000 K are independent on the occurrence of hydrogen shell 
flashes during the previous evolution. Hydrogen burning re- 
mains the main energy source along the co oling branch down 
to very low effective temperatures (see Sect. 4.3). Because He- 
WDs suffering from hydrogen flashes evolve back to the RGB 
regime one has to account for high mass-loss episodes due to 
Roche lobe overflow although the time spent away from the 
cooling branch is small compared to the cooling time itself. 
Our calculations show that the cooling is not affected by such 
high mass-loss episodes (see Sect. p~4l ). 



Fig. 3. Same as Fig. ||, but as a function of ipost-RGB- The age 
at the very left of the diagram corresponds to effective tem- 
peratures of about Toff ~ 15000 . . . 18000 K, the curves end 
at temperatures of T^s w 4000 . . . 5000 K. We adopted t = Q 
when the models pass Toff = 10000 K for M < 0.2 Mg and 
Tcs = 5000 K for M > 0.2 Mq after leaving the RGB. 



4.2. The main flashes 

Figures ^ and ^ show the complete evolutionary tracks for the 
sequences with M = 0.234 M© and 0.259 M©. While the evo- 



We did not find any unstable hydrogen-burning in the se- 
quences with masses M = 0.179, 0.195 and 0.414 Mg. Only 
for M = 0.234 and 0.259 Mq major hydrogen sheU flashes 
(^Hyd.max 10^ L©) developed with concomitant extended 
loops in the HRD. Only a temporal slight increase of the CNO- 
luminosity on the cooling branch was found in the sequences 
with M ^ 0.300 and 0.331 Mq. 

The restricted mass range for the occurrence of hydro- 
gen flashes agrees with earlier results, e.g. Kippenhahn et al. 
(1968), Gianonne et al. (1970) or Casteflani et al. (1994). 
We note that Webbink (1975) suggested a lower mass limit 
of M w 0.206 M©, in good agreement to our findings. The 
0.298 Mq model of Iben & Tututkov (1986) is fufly unsta- 
ble, our 0.300 Mq model only marginally. The calculations of 
Castellani et al. (1994) suggest that the upper boundary for the 
flash range is higher for lower metallicity. Therefore, their re- 
sult of M^^^^^^ < 0.33 Mq for Z = 0.01 is consistent with 
our calculations. Complementary calculations where equilib- 
rium rates for the pp-chains and the CNO-cycle were used in- 
stead of the nuclear network show three hydrogen shell flashes 
for M — 0.227 Mq and one strong hydrogen shell flash for 
M = 0.310 Mq. Because we found no hydrogen flashes in our 
standard sequences with M = 0.300 Mq and M = 0.195Mq 
we conclude that the mass range for the occurrence of instabili- 
ties depends on uncertainties due to different input physics (e.g. 
diffusion as in Iben & Tutukov (1986)), and that the boundaries 
M^J^^^. = 0.21 Mq and M^^^^ = 0.30 Mq may be uncertain 
by about 0.01 Mq. 

According to our calculations the number of hydrogen 
flashes depends on the mass: We found one major hydrogen 
flash for M ^ 0.25 Mq, two flashes for 0.25 ^ M/Mq Z 
0.23 and three flashes for 0.23 Z M/Mq Z 0.21. Thus, we 
conclude that unstable hydrogen burning in He-WDs caused 
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Fig. 4. HRD with a complete evolutionary track for the He- 
WD sequence with M — 0.234 Mq. The dashed line shows 
the 1 Mq track which was used to generate the post-RGB 
model. Two main flashes occurred, each followed by weaker 
subflashes (see also Fig. ^. 

lution of the M = 0.259 Mq model is characterized by only 
one major flash resulting in an extended loop in the HRD (see 
Fig. H), the model with M = 0.234 Mq experiences two strong 
flashes (see Fig. 0). 

After mass loss terminates the RGB evolution, the remnants 
evolve at almost constant surface luminosity to higher effective 
temperatures towards the white dwarf cooling branch. In this 
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velope mass and the thermal timescale of the burning shell, 
and ranges from Atonsct ~ 4 • 10'^ yr for M « 0.23 M© to 
Aionsot ~ 4-10'' yr for M « 0.30 Mq. For example, the model 
with M = 0.259 Mq gives Aionsct ~ 1-7 • lO'^yr, that with 



4.6 4.4 4.2 4 3.1 
logTeff[K] 
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Fig. 5. Same as Fig. |4|, but for the M = 0.259 M© sequence. 
Here only one major flash was found. Properties of the models 
labeled by capital letters are given in Tab. g. 

Table 2. Data for the marked points in Fig. |[ The age was 
set to zero at Toff = 10000 K of the post-RGB evolution. The 
different evolutionary stages A-G are discussed in the text. 

Model tpost-RGB/loVr log(reff/K) log(L/LQ) 



A 


1.469119 


4.2865 


-0.5726 


B 


2.803540 


4.3678 


-0.0754 


C 


2.803641 


3.9532 


-1.4070 


D 


2.803643 


3.8488 


-1.0786 


E 


2.803658 


3.8413 


-0.2007 


F 


2.803744 


4.3272 


2.0633 


G 


13.003598 


4.2859 


-0.8025 



0.234 Mq only Atonsot ~ 5 ■ 10'' yr. At the onset of the 

1.2-lO^yr 



M 

second flash for M = 0.234 M© one gets Afonsot 
due to the further reduced envelope mass. 

The increasing energy release during the flash development 
causes a steep temperature gradient in the vicinity of maximum 
energy production and the formation of a pulse-driven convec- 
tion zone wefl inside the hydrogen burning shell (beyond point 
B). The situation is displayed in Fig. |[ At t « 2.80357 • 10^ yr 



3.4 




Fig. 6. Development of the pulse-driven convection zone 
(intermediate grey shaded region) during the hydrogen-flash 
phase of the M — 0.259 M© sequence as a function of time. 
Point B from Fig. | is located at t « 2.80354 ■ lO'^ yr. Also 
shown are the lines of maximum energy generation due to hy- 
drogen burning, Cmax, the regions where e^ax has dropped to 
1% (dark grey shaded) and 0.1% (light grey shaded) respec- 
tively, and the degeneracy parameter ifj (right y-axis) at the 
point with e = Cniax- 



phase the total luminosity is almost entirely supplied by nuclear 
burning due to CNO cycling (L w Lqno)- The contribution 
due to contraction is negligible. The situation changes when the 
star reaches the cooling branch. Now, the temperature within 
the hydrogen burning shell becomes too low to support further 
CNO cycling. Lqno drops significantly, and contraction sets in 
until pp-burning takes over the main nuclear energy production. 
Still Lnyd 3> ig holds in this early cooling period (see e.g. Fig. 
[n] and Fig. 5 in Paper I). The phase of unstable burning starts 
at typical effective temperatures of T^s « 20000 . . . 25000 K 
on the cooling branch (close to point A in Fig. |]). At this 
stage of evolution the energy production due to CNO cycling 
rises again and finally exceeds the pp-contribution by large 
amounts. In Fig. ^ this part roughly coincides with the loop 
between point A and B in the HRD. The duration for this 
period of evolution (A^onsct) depends on the remaining en- 



the convective shell establishes right above the locus of max- 
imum energy production due to hydrogen burning. The de- 
generacy inside the burning shell is quite moderate as can 
be seen from the degeneracy parameter ip at maximum en- 
ergy generation (for the defintion of see e.g. Kippenhahn 
& Weigert 1990). During the onset of the unstable burning, 
V'«-0.5...-l. 

Once the convection zone is fully established, the evolution 
is rapidly accelerated (roughly between point B and D). The 
typical time scale is now of the order of a few decades. The lu- 
minosity contribution due to nuclear burning increases several 
orders of magnitude and the convection zone grows until it ex- 
tends to the stellar surface (see Fig. |, t « 2.80354 . . . 2.80364- 
10^ yr). Although Lnyd increases to about 10^ L© in this phase 
(near point C in Fig. ||), the surface luminosity drops by almost 
2 orders of magnitude because the increased energy production 
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is largely overcompensatedby the energy loss due to the expan- 
sion of the envelope. This expansion roughly doubles the stellar 
radius and leads to a complete lifting of degeneracy in the shell 
(sharp drop of For instance, model B has R « 0.05 R©, 
and beyond point C one finds i? « 0.1 Rq. 

The maximum hydrogen luminosity reached during the 
flash is supplied by pp-burning, although the onset of the insta- 
bility is triggered by CNO cycling (see below. Fig. For M = 
0.259 Mq, LHyd.max ~ 1-3 • IO^Lq. For M = O.234M0 
both flashes are of comparable strength with Lnyd.max ~ 
2.7 • 10^ L0 for the first and Lnyd.max ~ 2.3 • 10^ for the 
second flash. The time span between both flashes amounts to 
1.8- 10^ yr (see also Fig. 0). 

Beyond point C the most luminous part of the flash insta- 
bility has passed and Lnyd decreases while Lg increases. The 
surface luminosity is again enhanced due to the increase of 
Lg. At point D (Tos « 7000 K, see Fig. |]) the upper bound- 
ary of the pulse-driven convection zone reaches the surface (at 
t « 2.80364 • 10^ yr in Fig. |6|). Due to the convective trans- 
port of hydrogen from the surface layers down to the hydro- 
gen burning region the surface hydrogen abundance is reduced 
by about AX « 0.06 when surface convection establishes. 
The small blueward evolution in the evolutionary track close 
to point D is related to these composition changes in the en- 
velope (see Iben & Tutukov 1986). The compositional change 
is consistent with the results of the flash model of Kippenhahn 
et al. (1968) who found AX w 0.08. Iben & Tutukov (1986) 
found a larger amount of AX « 0.2, probably due to their 
consideration of gravitational and chemical diffusion leading 
to different chemical profiles. 

Between point D and E in Fig. |] the lower boundary of 
the pulse-driven convection moves upwards, and at point E the 
convection zone (and therefore surface convection) vanishes at 
all. Overall, the convection zone exists for « 1200 yrs, and for 
« 150 yrs it extends up to the surface. Beyond point E the 
contraction of the inner regions of the hydrogen shell resumes, 
while the surface layers react by expansion, resulting in a red- 
ward motion in the HRD bringing the He-WD almost back to 
the RGB region. Finally, contraction seizes the surface layers 
and the star evolves back to higher effective temperatures to- 
wards the cooling branch again. 

Around point F just before re-entering the cooling branch, 
the so-called subflashes develop. Here, contraction initiates an- 
other smafl increase of Lnyd to sa 30 . . . 200 L© (compared 
to « 10^ Lq during the major flash) leading to a temporary ex- 
pansion of the outer layers of the star and only minor circle-like 
excursions in the HRD. 

For illustration. Fig. || shows the evolution (arbitrary zero 
point) of the luminosity contribution due to hydrogen burning, 
LHyd, and the gravothermal luminosity, Lg, during the flash 
phase of the M = 0.234 Mq sequence. The large plot shows 
the evolution of the first subflash episode of the first hydro- 
gen flash phase and the small inlet that of the second phase. 
The peak at i « 0.4512 • 10^ yr which is not resolved within 
the plot range comes from the major flash, the second peak at 
t « 0.4518 • 10^ yr with inyd.max ~ 250 L© followed by a 



short period with Lg < 0, i.e. expansion, marks the subflash. 
This first subflash is then followed by another slight increase in 
Lnyd at t w 0.4528 • 10 '' yr with LHyd,max ~ 20 Lq causing 
the little circle in the track at log L/L© « 1.2 (see Fig. Eh. 



J" 




Age [in 10 yrs] 

Fig. 7. Evolution of surface and hydrogen luminosity, L and 
Lnyd, and the contributions due to CNO cycling and pp- 
burning as a function of cooling age for M = 0.234 M©. The 
figure shows the age range where the hydrogen shell flashes 
develop. Age t ~ is given by the first turning point after the 
track has entered the cooling branch (see point A in Fig. ^) and 
corresponds to a post-RGB age ipost-AGB ~ 4.62 • 10^ ■ 



yr. 
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0.454 
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Fig. 8. Evolution of inyd (dashed line) and Lg (solid line) 
as a function of time for the 0.234 M© sequence during flash 
phase: Large plot: First flash phase; inlet: second flash phase 
(same plot labels). For the definition of i = see Fig. ^ 

After the second major flash (again unresolved in the time 
range of this diagram at t w 2.2757 • 10^ yr) a local maxi- 
mum in ZvHyd with inyd.max ~ 30 L© and slight expansion is 
found at t « 2.2762 • 10'' yr (i.e. about 5 • 10^ yr after the main 
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flash) corresponding to the small loop in the evolutionary track 
at logL/L0 Rs 1.5. 

As previously mentioned the instability in He-WDs is 
caused by the fading CNO-luminosity on the cooling branch 
around T^s w 22000 . . . 25000 K. After pp-burning became 
the dominant contribution of hydrogen burning, Lcno again 
increases on a typical timescale of a few 10^ to 10^ yr Figure ^ 
shows the situation for the two hydrogen flashes occurring for 
M = 0.234 M0. Beyond t « 4.6 • lO'^ yr and < « 5.6 • 10^ yr 
LcNO^ rises again from Lcno ~ 0.1 to a few 10^ Lq. With on- 
going unstable burning the outer layers of the He-WD expand 
and slightly cool the shell so that Lcno drops again. The peaks 
in Lnyd (~ 10^ Lq) are due to pp-burning. 

To investigate the dependence of the mass on the occur- 
rence of hydrogen flashes we applied the criteria for unsta- 
ble burning discussed in Sect. 2A. As already seen in Fig. || 
the degeneracy in the hydrogen shell is quite moderate but far 
from being negligible. We used the criterion of Kippenhahn & 
Weigert (1990) to account for the possible effect of degener- 
acy. The criterion states that nuclear burning in a shell at radius 
i?shcn and of thickness AR becomes unstable if 

45 

1 - Vad • JT-^ > . (1) 



4a 



AR. 



Vad is the adiabatic temperature gradient and a and 5 are given 
by a = (f^Tp) ^ ^ ^ ( ainr ) ■ comparison, we 

also used the criterion of Schwarzschild & Harm (1965) which 



predicts unstable burning if 



AT 



> 



and 



AR 



<-,-\Q\ 



(2) 



Tshcll I' ^sh 

with Q « —8 ... — 4 and AT being the temperature differ- 
ence between lower und upper shell boundary. A main problem 
when dealing with these criteria is the definition of quantities 
as, for example, representative temperatures and pressures in 
the shell or the typical radial extent of the shell. 

However, although we find the instability criterions to be 
fulfilled in several models the strict application of the afore- 
mentioned criteria alone cannot explain the restricted mass 
range for the occurrence of hydrogen flashes. Some general 
aspects on the model properties important for the study of un- 
stable burning can be seen from Fig. ^ where different quanti- 
ties of layers at maximum energy generation are shown along 
the cooling branch evolution for the non-flash sequences with 
M = 0.195 and 0.414 M© and the flash sequence with M = 
0.259 M0 covering most of the He-WD mass range. 

As seen before for M — 0.234 Mq, degeneracy of the shell 
layers is moderate (lower left panel) in the upper part of the 
cooling branch evolution. For M ~ 0.259 Mq ,■(/;« —0.5 ... 
at the beginning of unstable burning (i « 2 • 10^ yr). For 
M — 0.414 Mq degeneracy is slightly lower {ip « — 1, 
t 2 ■ 10^ yr) at the corresponding evolutionary stage, and 
for M = 0.195 Mq sHghtly higher (ij; « . . . + 0.5, t w 
5 • 10* yr). For an ideal, non-relativistic gas one obtains for the 
non-degenerate limiting case a = 6 = 1, whereas strong de 
generacy {ip ^ 1) leads to a « |(1 + and 6 
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Fig. 9. Different quantities as a function of ipost-RGB for the 
sequences with AI — 0.195 and 0.414 (no flashes) and 
M = 0.259 M0 (one flash): Upper left: Temperature at the 
maximum of energy generation, e„iax, in the hydrogen burning 
shell. Upper right: Radial extent of the hydrogen burning shell. 
Lower left: Degeneracy parameter t/j at e = Emax- Lower right: 
Temperature derivative of the energy generation rate, ex, at e = 



with iio and /ic being the mean molecular weights per ion and 
electron, resp., (see Kippenhahn & Weigert 1990) |. Thus, in- 
specting Eq. (|l]), from ^ alone we would conclude that, if at 
all, degeneracy favours hydrogen shell flashes for lower He- 
WD masses. 

The temperature sensitivity of nuclear burning along 
the cooling branch can be seen from the quantity ex = 
d\ne/d\nT{= viore oc T") which is shown in the lower 
right panel of Fig. ^ When the flash develops for M = 
0.259 Mq, one finds ex ~ 12.5. For the other masses we have 
ex w 15 (M = 0.414 Mq) and ex « 11 (M = 0.195 M©). Be- 
cause unstable nuclear burning is favoured for higher temper- 
ature sensitivity of the burning (see Eq. (||)) the occurrence of 
flashes is favoured for heavier white dwarfs. Thus, combining 
the effects of degeneracy and large temperature exponents in 
the energy generation rate might, in principle, explain why high 
and low-mass He-WDs do not suffer from hydrogen flashes. 



3 Ale J_ 

2 /io V' 



Note that radiation pressure is negligible and Vad does not de- 
pend on in the non-relativistic regime. 
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However, it is most likely the radial thickness of the shell 
(see upper right panel in Fig. ^ which has the most important 
effect on the occurrence of hydrogen flashes. At the onset of un- 
stable burning we found AR « 10"^ R© for M = 0.259 Mq 
(the shell's borders are taken at the point with e = 10^'^ • emax)- 
For M = 0.195 Mq, AR is approximately a factor of two 
larger and for M = 0.414 Mq a factor of two smaller. For low- 
mass He-WDs degeneracy of the shell is of minor importance. 
Also, the thickness of the shell is too large as to allow hydro- 
gen shell flashes. On the other hand, for high-mass He-WDs 
the shell is so thin that the thermal cooling time is smaller than 
the typical time scale for the onset of the instability. Hence, for 
He-WDs with M 0.3 M© an instability might be initiated 
but the fast cooling of the shell prevents a hydrogen flash (see 
also next section). 



4.3. Final cooling 

After the He-WD models with M = 0.234 M© and 0.259 M© 
have completed their major hydrogen flashes a final onset of 
a flash results in a local hook in the cooling track (see Figs. Q 
and ^ point G). Again this increase of luminosity due to hy- 
drogen burning is caused by CNO cycling. This can be seen 
from Figs. |l^ and |ll] where the evolution of different lumi- 
nosity contributions for the two flash sequences are plotted as 
a function of tpost-RGB for the whole post-RGB evolution. 



O 




10=^ 
Age [in yrs] 



Fig. 10. Evolution of different luminosity contributions (see 
plot labels) as a function of tpost-RGB for the sequence with 
M = 0.234 M©: Lpp! luminosity due to pp-burning; Lcno- 
luminosity due to CNO cycling; LB_yd' complete luminosity 
due to hydrogen burning; L: surface luminosity; L^: contribu- 
tion of gravothermal luminosity due to contraction {Lg > 0); 
Ln luminosity contribution due to neutrino losses, t — is 
adopted as in Fig. I (see also Fig. | for i = 3 . . . 8 • lO^yr). 



LcNO shows a local maximum at tpost-RGB ~ 1-23 • 10 yr 
for AI — 0.234 M©, approximately 5 • 10^ yr after the last ma- 
jor flash. The corresponding point for M = 0.259 M© is at an 
age of ipost~RGB ~ 1-40 • 10^ yr, which is about 10* yr after 
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Fig. 11. Same as Fig. |l^ but for the M = 0.259 M© sequence. 



the major flash. At this point of evolution hydrogen shell burn- 
ing is already dominated by pp-burning (L 



pp 



20 • LcNo). 

but the difference between both contributions is rather small 
compared to the later evolution (see Figs. 10 and [iT] for t ^ 1 
Gyr). Because the characteristic cooling time is now smaller 
than the typical timescale for the onset of unstable burning an- 
other strong hydrogen flash is prevented. Such onsets of un- 
stable CNO-burning have also been found in the sequences 
with 0.300 M© and 0.331 M©, resulting in similar hooks on 
the cooling track (see Fig. 1 in Paper I). 

The further evolution in both flash sequences is comparable 
to the one without hydrogen flashes (see for instance Fig. 4 in 
Paper I). It is characterized by quiescent hydrogen shell burn- 
ing via the pp-chains until an age of about 10 Gyr. Then hydro- 
gen burning terminates and contraction determines the surface 
luminosity evolution and final cooling (see Figs. |l^ and [TTI ). 

Finally, we like to note that the mass-radius-relation of He- 
WDs is only affected by hydrogen shell flashes for Toff « 
20000 . . . 25000 K: The high burning rates associated with the 
flashes lead to a higher envelope consumption compared to se- 
quences without flashes, resulting in slightly smaller radii and 
enhanced evolutionary speeds on the upper part of the cooling 
branch. Besides, since the timescale of the instabilities is short 
compared to the characteristic cooling time (a few 10^ yr com- 
pared to Gyr) and most of the time during a flash is spent close 
to the cooling branch, the influence of hydrogen shell flashes on 
the mass-radius-relation is restricted to the onset phase where 
the changes in radius are moderate. 

4.4. Roche lobe overflow during the flashes 

During the hydrogen shell flashes the expansion of the outer- 
most layers (from a few 0.01 solar radii to a few solar radii) 
forces the He-WDs to evolve back into the RGB domain in 
the HRD (see Figs. ^ and Due to the binary nature of He- 
WDs, in principle, mass transfer to a companion (in most cases 
another white dwarf or a pulsar) should be taken into account 
in this particular evolutionary phase. This has been considered. 
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for instance, in the calculations of Iben & Tutukov (1986) lead- 
ing to a significant removal of envelope matter and a faster ex- 
haustion of hydrogen burning. 

Besides the sequences already discussed we thus also in- 
vestigated the evolution of the flash model sequence with M = 
0.259 M0 considering high mass-loss episodes due to Roche 
lobe overflow. Because we were only interested in the general 
influence on the final cooUng properties of our models we used 
a rather simple algorithm to account for Roche lobe overflow 
from the He-WD to a companion: For a given Roche radius 
^Rochc (in our case 2, 3 and 5 R©) and actual model radius 
R, we increased 77 in the Reimers mass loss formula by a factor 
of 2 • R/Rjxochc for R > i? Roche- Otherwise, rj was set to its 
standard value of 0.5 used for all post-RGB calculations. 

Figure |l^ presents the evolutionary tracks in the HRD 
during the flash for M = 0.259 Mq and different assump- 
tions of the Roche radius. The models quickly reach the line 
with R = i?Fiociio- For comparison, also the track without 
Roche-lobe overflow is shown. During the adjustment phase, 
after the models have passed R = i?Rochc, mass-loss rapidly 
rises up to A/high ~ lO~^M0/yr and forces the models to 
evolve back to i? w i?Roche- Finally, M « 5 • lO^'^Mg/yr 
when the models evolve along at i? « i?Rochc- When Roche 
lobe overflow ends (R < i?Rochc) mass loss has dropped to 

M = A'/Reimcrs « 10"" Mg/yi . 

After the end of Roche lobe overflow and the subflash phase 
the tracks rapidly merge with the one calculated without Roche 
lobe overflow at the beginning of final cooling evolution. As 
Fig. |l^ shows the subflash excursion is the less pronounced the 
smaller -RrocIio- 

The convergence of the tracks shows that the cooling prop- 
erties of He-WD models are not seriously affected by the Roche 
lobe events because similar tracks imply similar radii and thus 
similar envelope masses which determine the cooling evolu- 
tion due to hydrogen burning. This is confirmed by Figs. |l3|, |l4| 
and |l5[ In Fig. |l3|the evolution of the envelope mass, Mem, is 
given as a function of Toff for the tracks from Fig. |l^ in Fig. |l5] 
the corresponding evolution of surface luminosity as a function 
of ipost-RGB is plotted. Figure |l^ shows the reduction of en- 
velope mass for all sequences as a function of ^post-RGB close 
to R = i?Rochc and beyond. Figures [l^ and |l^ indicates that 
~ 6 • 10^* Mq of the envelope mass is lost in the non-Roche 
model during the flash phase (taken from the first local maxi- 
mum in Tcff to the last one). Note that the temporal increase of 
Afonv in Fig. 13 comes from the pulse driven convection zone 
by mixing hydrogen-rich material somewhat below the layers 
with X = 0.35 which determine the border between core and 
envelope (cf. Fig. 

A significant reduction of Mcnv sets in when dominant hy- 
drogen burning has again established at the end of the sub- 
flash phase ( see Fig. |l3] at Toff « 27000 K and Fig. [ij for 
2.81 ^ t/lO^yr ^ 2.95). This envelope reduction is slowed 
down when the model enters the cooling branch (see Fig. 
at Teff,max ~ 33000 K and Fig. [l^ for t ^ 2.95 • 10^ yr). 
The envelope is eroded by mass loss and shell burning. This 
can be written as Afcnv = A/wind + Afcoro with the mass loss 



rate A/wind and the core growth rate Mcorc- The core growth 
rate is given by A/corc — x^e^ ^''^ hydrogen content X 
and the gain of energy due to hydrogen burning per mass unit, 
«6.3-10iScrgg-i. 

For sequences without Roche-lobe overflow Afonv is 
mainly determined by shell burning, i. e. A/cnv ~ Afcore- 
With a mean value of Lnyd ~ 30 L0 and X « 0.7 it is 
A/coro « 4.3 • 10-1'' M© yr-i. With At « 1.3-106 yr (see Fig. 
|l^ one obtains AAfonv ~ Mcoxc ■ At « 5.6 • 10""^ Mq in good 
agreement with the value obtained from Figs. |l3| and |T^. 

For sequences with Roche-lobe overflow the situation is 
different. Here, the consumption of envelope mass occurs pref- 
erentielly due to the high mass-loss phase when the models are 
close to their Roche limit (Teg « 8000 . . . 15000 K, see Fig. 
pj[ ). For instance, for the sequence with i?Rochc — 3 Rq it is 
AAfonv ~ 4.5-10^'^ M0, corresponding to a mean mass loss of 
A^wind ~ lO"^M0yr~i lasting for At « 40 yr. Afterwai'ds 
Afcnv is further reduced by AAfcnv ~ 2 ■ 10^^ Mq due to nu- 
clear burning (see Fig. n^. 



S 

'o 




25 20 15 

T^ff [in 1000 K] 

Fig. 13. Evolution of Afcnv as a function of Tofi for tracks with 
Af = O.259M0 and different Roche radii (see labels). The 
smaller f?Rochc the faster Afcnv is reduced, but when the star 
is back on the cooling branch the residual envelope mass is al- 
most the same for all Roche tracks as it is for the non-Roche 
lobe track. 

In the end, the envelope masses are almost the same when 
the different models enter the final cooling branch (Afcnv ~ 
2.05 - 10^^ M0, see Fig. 14), i.e. the residual envelope mass 
and the corresponding cooling evolution are rather independent 
from the evolutionary history. 

On one hand, the latter result is consistent with the results 
of Iben & Tutukov (1986). In the case of Roche lobe over- 
flow the high mass loss determines the evolutionary speed, 
whereas otherwise the evolution is characterized by the thermal 
timescale of the envelope. On the other hand, in their 0.298 Mq 
model the high mass loss causes a considerable loss of enve- 
lope mass (« 1.5 - lO^'^M©) and prevents hydrogen burn- 
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Fig. 12. HRD with tracks for A/ = 0.259Mq and different Roche radii (2, 3 and 5 R©). For R > i?Rochc mass loss was adjusted 
as to provide R « -RRochc as long as the star does not evolve back to the white dwarf domain. Lines of constant radius are also 
plotted. The arrows indicate the course of evolution. The inlet in the lower right corner is a blow-up of the part close to the Roche 
radii. 



ing from becoming ever dominant again. Consequently, their 
model cools down quite rapidly. The cooling time is of the or- 
der of a few 10® yr whereas our models are controlled by resid- 
ual hydrogen burning leading to cooling times of the order of 
Gyr (see Figs. and [ill). 

Thus, from our calculations it appears that different mass 
loss histories during the phase of the hydrogen shell flashes 
do not influence the final cooling of He-WDs significantly. 
This result seems to be consistent with the calculations of Kip- 
penhahn et al. (1968): Their M = 0.264 M© model evolves 
through one major hydrogen shell flash and one strong sub- 
flash, both bringing the white dwarf radius above the critical 
Roche radius (points Q to R and T to U in their Fig. 1), result- 
ing in a total mass loss of about lO^'^ M© which is comparable 
to the 6 ■ lO^** Mq of the present calculation with only one high 
mass-loss episode. After a rapid evolution back to the cooling 
branch the evolution is slowed down, and at log(L/L0) « —2 
(point W in their calculation) a cooling age of approximately 
2.55 Gyr is reached, comparable with our M — 0.259 M© 
model which has at the same position tpost-RGB ~ 1-85 Gyr. 
On the other hand, by assuming a linear correlation between 
log L /Lq and log tcooi as a first approximation, the cooling age 



of the model of Iben & Tutukov (1986) can be estimated to be 
onlyicooi ~ 2.3-10® yr at log(L/L0) « —2 (using their points 
S and T, see their Fig. 1 and Table 1). This value is almost one 
order of magnitude below the one given by models which allow 
hydrogen burning to continue. 

5. Summary 

We have calculated evolutionary models of low-mass white 
dwarfs with helium cores to study in detail their cooling be- 
haviour and their evolution including phases of thermally un- 
stable hydrogen burning (hydrogen shell flashes). We found 
that the occurrence of these thermal instabilities is restricted 
to the mass range 0.21 ^ M/Mq ^ 0.3. It is noteworthy that 
a sufficient temporal resolution is essential to avoid numerical 
fluctuations in the energy output of the shell during early cool- 
ing branch evolution. 

The flashes occur during the fast cooling of the shell at the 
beginning of the cooling branch evolution. This especially af- 
fects the lower part of the shell region where CNO cycling is 
the dominant contribution to hydrogen burning. CNO-burning 
is temperature sensitive enough to cause thermal instabilities 
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Fig. 14. Evolution of envelope mass as a function of ^post-RGB 
for tracks with M = O.259M0 and different Roche radii (see 
labels). The plot shows the reduction of A/onv during and im- 
mediately after Roche-lobe overflow. The inlet shows the evo- 
lution of il/cnv on a larger scale. 




10^ 10^ 10^° 

Age [yr] 

Fig. 15. Evolution of surface luminosity i as a function 
of ipost-RGB for tracks with M ~ O.259M0 and different 
Roche radii (see labels). The models suffering a high mass-loss 
episode show the same cooling behaviour as the non-Roche se- 
quence. The inlet shows the evolution of L during the subflash 
phase. 



This result also holds when Roche-lobe overflow of the He- 
WD to a companion is considered during the expansion into 
the RGB regime. One of the main features related to hydrogen 
shell flashes is the change of envelope composition due to a 
pulse-driven convection zone which temporarily extends to the 
outermost layers of the star. The resulting drop in hydrogen is 
AX « 0.06. 

While the typical duration of a hydrogen flash is short com- 
pared to the cooling time of He-WDs (10^ yr compared to Gyr) 
and He-WDs evolve close to their cooling tracks for most of 
the time during a flash, the influence of these unstable phases 
on the mass-radius-relation of helium white dwarfs is moderate 
and restricted to effective temperatures of 20000 . . . 25000 K. 

Acknowledgements. F.H. acknowledges funding by the Deutsche 
Forschungsgemeinschaft (grant La 587/16). 
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if two conditions are fulfilled: The shell is thin enough (this 
is obviously not the case for M/M0 ^ 0.21) and the cooling 
time is large enough as to avoid an extinction of the shell be- 
fore the instability is fully established (this does not hold for 
Af/M0 ^ 0.3). The final cooling is not affected from possi- 
ble flash events, i.e. hydrogen shell burning establishes again 
as dominant energy source at the end of the flash episode as in 
the case of non-flash models. 



Appendix A: Remarks on the appropriate temporal 
resolution 

The question of sufficient resolution always arises when one 
deals with instabilities in evolutionary models. As mentioned in 
Sect. 2.2 thermal instabilities in He-WD model calculations can 
obviously be missed if too large time steps are used. Large time 
steps may also cause fluctuations in the shell energy output. 
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Fig.A.l. Hertzsprung-Russell diagram for an evolutionary 
track with 1\I = 0.195 M0 and different temporal resolutions 
(see labels). The upper left panel shows the final track with 
the appropriate resolution giving a smooth course of the hydro- 
gen luminosity contribution Lnyd (see Fig. A. 2). For the other 
tracks deviations caused by local fluctuations in Lnyd (see Fig. 
A. 2) increase with increasing maximum allowed time step. 



pled the evolutionary time step to changes in Lnyd by reducing 
Aimax by a factor of 2 if inyd changes by more than 5 %. 
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Fig. A.2. Evolution of the luminosity contribution due to hy- 
drogen shell burning (Lnyd) for M = 0.195 Mq with different 
temporal resolution (see also Fig. A. 1 ). 



The four panels in Fig. A. 1 show a part of the evolutionary 
track of our He- WD model with M = 0.195 Mq in the vicinity 
of the turn-around point, Toff = Toff .max with different tempo- 
ral resolution. The corresponding evolution of the luminosity 
contribution due to hy drog en shell burning, Lnyd, is shown in 



the four panels of Fig. A. 2 



The track with sufficient temporal resolution (upper left 

lO^'yr) shows no perturba- 



panel in Fig. A.l, At^ax = 2 
tions due to fluctuations in the shell energy production (smooth 
curve in Figs. A.l and A. 2). When the maximum time step is 
increased to Aimax — 5 • lO'^yr small perturbations occur in 
Lnyd and in the ti-ack (for -0.4 < log L/LQ < -0.25, see 
upper right panel in Fig. AA ). When the maximum time st ep is 



further increased by a factor of 2 (lower left panel in Figs. A. 1 



and |A.2| ) the fluctuations in Lnyd and in the track are much 
more prominent but the evolution stabilizes after a final hook 
in the track at logL/L0 —0.5. In this case the shape of the 
tracks corresponds with those of e.g. Webbink (1975) or Sarna 
et al. (1998). Finally, for Aimax = 5 • lO^yr the calculations 
became numerically unstable. 

Therefore, we selected Aimax ^ 2 • lO^yr as the maximum 
time step for our calculations. Additionally, to handle the large 
local changes in the luminosity budget of the models we cou- 



